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0m . Abstract 

CU ■ 

I , We consider the stochastic wave equation with multiplicative noise, which 

i-Q ■ is fractional in time with index H > 1/2, and has a homogeneous spatial 

covariance structure given by the Riesz kernel of order a. The solution 
is interpreted using the Skorohod integral. We show that the sufficient 
condition for the existence of the solution is a > d— 2, which coincides with 
the condition obtained in [4|, when the noise is white in time. Under this 
condition, we obtain estimates for the p-th moments of the solution, we 
^ ' deduce its Holder continuity, and we show that the solution is Malliavin 

l/~) , differentiable of any order. When d < 2, we prove that the first-order 

Malliavin derivative of the solution satisfies a certain integral equation. 
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1 Introduction 

In the present article, we consider the following Cauchy problem: 



— y(*,x) = Au(t,x)+uW(t,x), t>0,xeR d (1) 

u(0,x) = u Q (x), x e K d 
^(0,x) = «i(a;), xeR d , 

where uq and u\ are deterministic functions, and W is a zero-mean Gaussian 
process which is fractional in time, with Hurst index H > 1/2, and homogeneous 
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in space, with covariance kernel / (to be specified below). In other words, 
W = {W(<p); <p £ HV} is an isonormal Gaussian process, defined on a complete 
probability space (fi,.F, P), with covariance: E[W(<p) W{ip)] = (<p,ij))wp. 

Throughout this article, T-LV denotes a Hilbert space (which may contain 
distributions in S'(R d+1 )), defined as the closure of the set £ of linear combi- 
nations of elementary functions l[o,t]xA,t >0,4€ i3f,(R d ), with respect to the 
inner product: 

(l[o,t]xA>l[o,.]xB)w ~Rn(t,s) / f(x-y)dxdy. (2) 

J A J B 

(Here, Bb{M. d ) denotes the class of bounded Borel sets in R d .) 

The notation appearing in ^ needs some explanation. Rn(t, s) denotes the 
covariance of the fractional Brownian motion of index if, and since we assume 
that H > 1/2, we have: 



R H {t 7 s) =0!H 

Jo Jo 



t r s 

2H-2 



\u — v\ 



dudv, 



where an — H(2H — 1). On the other hand, / : W l — >• M + is the Fourier 
transform in iS'(R d ), of a tempered measure /i on R d , i.e. 

f(x)4>(x)dx = f r<KZ)ii(dO, v^ e s(R d ), 



where F (/)(£) = J„ d e % ^' x 4>(x)dx is the Fourier transform of 0, and 5(R d ) is the 
space of rapidly decreasing C°°-functions on R d . 

In the present article, a solution of ([T]) is an adapted square-integrable pro- 
cess u = {u(t, x); t > 0, x £ R d } which satisfies the following integral equation: 

u(t,x) = w(t,x)+ / G(t-s,x-y)u(s,y)W(Ss,Sy), (3) 



where G is the fundamental solution of the wave operator, the stochastic integral 
is interpreted in the Skorohod sense, and w is the solution of the equation 
w t t — Aw, with initial conditions w(0, ■) = Uo, wt{0, ■) = u±. 

We are interested in the case uq = 1 and u\ = 0, and hence w(t, x) = 1 
for all t > 0, x £ R d . (In the case H — 1/2, this corresponds to the equation 
Uu — Aw + (u + 1)W with zero initial conditions.) 

The case of the heat equation with multiplicative fractional noise was treated 
in [S] , [B] and [T] , while the recent article [7] gives a Feynman-Kac representation 
for the solution. In the case of the heat equation, the key estimate which leads 
to a sufficient condition for the existence of the solution is: 

f e~ t ^/' 1 \^-r]\- a d^<C a4 t-^- 0i y 2 , foranyt>0,?7eR d , (4) 

Jv. d 

for any < a < d (see Lemma 6.1 of [5], or Lemma 3.3 of p]). 



The following estimate lies at the origin of our developments, being the 
analogue of ((4]), which is needed for the wave equation: 



: "' m Jr% - v\- a <% < G aA t a - d +\ for any t > 0, r, e 



* 2 ( 

Jm* ici 2 

for any d — 2 < a < d (see Lemma 13.11 below) . 

As far as we know, the only other study of the wave equation, in which the 
fractional noise enters the equation in a multiplicative way is [12] . which treats 
the case d = 1. The authors of [12] use a pathwise integral for interpreting 
the solution, instead of the Skorohod integral used in the present article, which 
makes it difficult to compare the results. 

The study of the wave and heat equations driven by a Gaussian noise which is 
white in time and homogeneous in space was considered by many authors, using 
martingale methods (see @], [5], Q3] , [3], and the references therein). These 
methods work for more general equations (in which the factor u multiplying 
the noise may be replaced by cr(u), for a Lipschitz function a), but fail in the 
fractional case. The method that we use in the present article is specific to the 
case <j{u) — u, in which the solution has a Wiener chaos decomposition whose 
kernels can be written down in closed form. 

The fact that the solution has a known Wiener chaos decomposition has sev- 
eral implications. More precisely, once we have a uniform bound for the moments 
of order 2 of the solution, we can pass to the moments of order p > 2, using the 
hypercontractivity of the Ornstein-Uhlenbeck semigroup. In this manner, we 
obtain estimates for the p-th moments of the increments of the solution (which 
yield its Holder continuity), and we show the Malliavin differentiability of the 
solution. Finally, assuming that d < 2, we prove that the first-order Malliavin 
derivative of the solution satisfies an integral equation, using a Hilbert-space- 
valued Skorohod integral. These results are valid for the heat equation as well. 

The present article leaves some open problems: 

1. For the existence of the solution, we could treat only the case of the Riesz 
kernel f(x) = c a .d\x\~ i - d ~ a ' > with < a < d, but other particular cases (e.g. 
f( x ) — Yli=i a Hi\%i\ *~ 2 with 1/2 < Hi < 1) would be similar. However, we 
do not know what is the sufficient condition for the existence of the solution, 
in the case of a general kernel /. Even in the case of the Riesz kernel, it is not 
clear if the sufficient condition a > d — 2 is necessary for the existence of the 
solution. These questions are open in the case of the heat equation as well. 

3. In the case of the white noise in time (see e.g. [11] V the argument which 
leads to the existence of the density of u(t, x) relies on the fact that Du(t, x) 
satisfies a certain integral equation, and its "HP-norm coincides with the norm in 
L 2 ((0, t);H), for a certain Hilbert space H of distributions in S'(R d ). Therefore, 
in the calculation of ||£>u(t, x)||-wp, one can bound from below the integral over 
(0, t) with the integral over a small interval (t, t — S). This bound is not justified 
in the fractional case. As we could not find an alternative argument to prove 
that ||Du(i, x)||w > a.s. (which would allow us to apply the Hirsch-Bouleau 
criterion), we could not show that the law of u(t, x) has a density, even if d < 2. 



2 Preliminary Results 

Intuitively, the solution of ([1]) should be given by a series of iterated integrals: 
u(t,x) = l + / G(t-ti,x-Xi)W(dti,dxi)+ 

/ G(t-t 2 ,x-x 2 )( / G{t 2 -h,X2-Xi)W(dh,dxi))W(dte,dx 2 )+.. 

JR d \J0 JR 2 J 

Since in dimension d > 3, Git, ■) is a distribution in R d , the product 

In \p\t X\ , . . . , &7i ; X n , L, X ) . — Lj[t t n , X X n JLr^6 n t n _i , X n X n — 1 J • • • 

G(i 2 - <i,*2 ~ Ki)l{o<ti<...<t„<t}> ( 5 ) 

has to be defined as the product of distributions, and one has to be careful with 
the definition of the iterated integrals above. 

The goal of this section is to take care of this difficulty, by tackling the 
following three problems: 

• (Subsection 12. ip For any t > 0, x £ R d and for any < t\ < . . . < t n < t, 

we give a meaning to f n (ti, •, . . . ,t n , -,t,x) as a distribution in S'(M. nd ), 
and calculate its Fourier transform. 



• (Subsection I2.2[) For any t > 0,x £ M d , we give a general criterion which 
ensures that f„(-,t,x) £ HV® n . 

• (Subsection ESJ Suppose that for any t > 0, x £ R d , /„(•,£, «) £ UV® n , 
and the series u(t, x) := l + J2 n >i In{fn{-,t,x)) converges in L 2 (fi), where 
I n denotes the multiple Wiener integral with respect to W. We show that 
u = {u(t,x);t > 0,x £ M. d } is a solution of (JlJ, in a sense which will be 
described below. 

2.1 The definition of the kernels f n (-,t,x) 

In this subsection, we give a rigorous meaning to the kernels f n (-,t,x). We let 
C^°(M. d ) be the space of infinitely differentiable functions on M. d with compact 
support, and V(R d ) be the space of (Schwartz) distributions on R d . 

Assume first that n — 2 and let < t\ < t 2 < t be arbitrary. We proceed 
to the formal calculation of the action of fiiti, ~,t 2 , -,t,x) on a test function 
4> = 0i ® 2 with 0i, 2 € C* °°(R d ): 



(f2(ti,-,h,-,t, x),4>) = / G(t-t 2 ,x-x 2 )(f> 2 (x 2 ) G(t 2 -ti,x 2 -xi)(f>i(xi)dxidx 2 

JTSL d 

G(t- t 2 ,x - x 2 ) 02 (x 2 )ipi(t 2 - h,x 2 )dx 2 



G(t -t 2 ,x — x 2 )tp 2 (t 2 ~ ti,x 2 )dx 2 
Viit 2 — ti,t — t 2 ,x) 



where 

V'i(-) = <M0> and ^2(s,-) = <M-Vi( s > - ) 

lf!(s, •) = ^l(-) * G{s, •), V2(Sl, S2, •) = ip2(si,-) * G(s 2 , ■), 

and * denotes the convolution with respect to the space variable. Similar formal 
calculations can be done for any n. 

Based on these calculations, for any < ti < ... < t n < t fixed, we let 
f n (ti, ■,..., t n , •, t, x) be the element of V(W nd ) whose action on a test function 
<f> = fa ® . . . ® fa with 4 e Cg°(R d ), is given by: 

(/n(*l, ■,■■ -,tn, -,t,x),4>) ■— ip n (t 2 -*1,*3 ~*2, •••,*- *»>#)> ( 6 ) 

where the pairs (^ , <£& ) are defined recursively for k = 1, . . . , n by the following 
relations: 

^fe(si,...,Sfc_i,-) = fe (-)^fc_i(si,...,Sfe_i,-) (7) 

( y 5 fc (si,...,s fc ,-) = ipk(si,.-.,s k -i,-)*G(s k ,-). (8) 

Note that Vfc(si,---,s fc -i,-) € C °°(]R d ) and ¥>fc(*i,...,*fc,0 G 5(R d ), since 
G(s,-) is a distribution with rapid decrease in M. d (see p. 245 of 14). The 
previous definition is extended to = fa ® . . . ® fa with 0,; G S(M. d ). 

The next result shows that the Fourier transform of f n (t\, •, . . . , t n , •, t, x) is 
a function in W ld , given by: 



Tf n (h ,-,... , t n , -, t,x)(£ u ..., tin) = e - ^+- + ^-*TG(t 2 - t u Qfo) 

■TO(*3 - *2, 0(6 + 6) ■ • ■ FG(t - t„, 0(6 + • ■ ■ + 60- (9) 

Proposition 2.1 For any < £i < . . . < £„ < t and for any h — hi ® . . . ® /i n 
mf/i /i, e C£°(R d ), we /law: 

(/„(<!, •,..., t„, •, t, x), h) = [ hfa, • . • , Z n )e- i & + " + M- x rG(t 2 -t 1 ,-)(Z 1 ) 



FG{U - fc, 0(6 + 6) • • -^G(t - i n , 0(6 + • • • + 6)^6 • ■ • #»• 

Proof: Note that := J"/i = fa <g> . . . ® 0„, where 4 := J"/i 4 e 5(R d ). By the 
definition of the Fourier transform in <S'(R d ) and (JSJ, we have: 

(Ffnih, •,..., t n ,-,t,x),h) = (/„(ii, •,..., t n ,-,t,x),fa) 

= (p n (t2-h,t 3 -t 2 ,...,t-t n ,x), (10) 

where (V>fc, <y9fc), 1 < fc < n are defined recursively by (O-©. 

We proceed to the evaluation of (f n (t 2 — t±, £3 — t 2 , ■ ■ ■ , t — t n , x). 

Step 1. For any si,...,s n _i G [0,i], we define recursively the following 
functions: g\ = hi, 



gk(si, ■ ■ ■ , Sfe_i,0 = /»fc * (fffc-i(si,. ■ • , Sk-2, ■)J 7 G(sk-i,-)), 2 < k <n. 



By induction on fc, 2 < k < n, it follows that gk(si, ■ • ■ , Sfc-i, •) € Cg 
(since J 7 G(s, •) is a G°°-function on R d ), and 



fffe(si, . . . , Sfe-1,%) = / h 1 (r] 1 )h 2 (r]2-Vi) ■ ■ ■ ^fe(%-%-i)-? r G ! (si, -)(?7i) • • • 

(fc-l)d 



J"G(s fc _i,-)(r?fc_i)dr?i ...dr\ k -\- (11) 

.Step D. We prove by induction on fc, 1 < A: < n, that: 



¥>fc(si, . ■ • , Sfe, •) = J"[5fc(si, ■ • ■ , Sfe_i,-)J"G(sfc, •)]■ (12) 

Before proving (O, note that: 

Fg*G( S ,-) = F(gFG(s,-)), Wg € G °°(R d ),V,s > 0, (13) 



since (J" 5 ) * G(s, •) = T\T~ x (Tg * G(s, •))] = T{gFG{s, ■)). 
For fc = 1, we have g\ = hi = T~ x 4>\ and 



<pi(s, •) = 0i * G(s, •) = J"g! * G(s, •) = T(giTG(s, •)), 

where we used (fT3"|) for the last equality. This proves (fT2")) for fc = 1. 
Suppose that (O holds for fc - 1. Then 

■tpk(8i,...,8k-i,x) = fe (a;)( y 3 fe _i(si,...,Sfc_i,x) 



= J r ft, fe (a:)J r [5fc_i(si,. ..,s fc _2,-)- F G(s fc -i,-)](a:) 
= T[hk * (fffc-i(si,. ■ • ,Sfc_2, ■)-7 r G(sfe_i,-))](a;) 
= J r .gfc(si,---,s fc -i,-)(a;) 

and 

^fc(si,...,s fc ,-) = ipk(si,---,s k -i,-) *G(sk, •) = ■? r S'fe(si,---,Sfc_i,-) * G(sfe, 

= ^"^(si, ■ ■ ■ , s fc _i, •)J"G(s fc , •)], 



where we used p3|) for the last equality. This proves (|12|) . 

Siep 5. Using (IT21 and (j2.ll) . we obtain that for any 1 < fc < n, 



(pk(si,...,Sk,x) = / e l,,k ' x g k (si,. . . ,s k -i,r]k)J r G(sk,-)(rik)dr)k 



e M7fe ' x /ii(77i)/i 2 (r? 2 - »7i) ■ • ■ M% - %_i).FG(si, -)M 



. .. J"G(s fe _i,-)(r7fc_i) J r G(sk,-)(Vk)dr)i ...dr]k 

where for the last equality we used the change of variables £i = 771, £j = r/j—rij-i, 
for 2 < j < fc. We now use ([TO]) . The conclusion follows using (Tl4l) with k = n, 
s n — t — t n and Si = i;+i — U for 1 < i < n — 1. □ 



2.2 The space UV® n 

In this subsection, we give a criterion for checking that an element ip £ r H,V® n , 
which can be viewed as a multi-dimensional analogue of Theorem 2.1 of [2]. 
This criterion is then applied to the case if — /„(-, t, x). 

For this purpose, for any T > fixed, wc define the multi-dimensional trans- 
fer operator K* H n by: 

n 

( K H,n 1 [o,t 1 ]x...[o,t„]){si, ■ ■ ■ , s n ) := JJ(^l[o,t i ])(s»), si,...,s n £ (0,T). 

f=l 

Let £c(0, T) be the set of all complex linear combinations of indicator func- 
tions l[o.t],i £ [0,T], and Hc{Q,T) be the closure of £c(0,T) with respect to 
the inner product: 

T r T 

2H-2, 



JO Jo 

The operator K* H n is an isometry between £c(0, T)® n and L 2 ((0, T) n ), which 
can be extended to l-Lc{Q,T)® n . In terms of fractional integrals, we have: 

(K* H ^)( S ) = (c^rr(F-i/2)"[ s ] 1 /2-H / ^-y2 ([u] H-i/ 20(u))(s)) 

where c^ = ( ^ (g _ 1 /g 2 _ 2g ) ) 1/2 > s = ( s i>- • • > s »)> [s] = Si • • . s n , and 

W,, ^(-( IU - ,rl/(n)dU 



is a multi-dimensional fractional integral of / e L 1 ((0,T)™), of order a £ (0, 1). 
For any function cj> £ Hc(0,T)® n , we have: 



a n H [ 0(u)0(v)[u-v] 2H - 2 dudv = 

J(0,T) 2 " 

^ / l^» 2 ({u] H -V 2 cj>(xi))( S )\ 2 \ H , n (ds), (15) 

J(0,T)" 

where d H = (c^) 2 r(H - 1/2) 2 and X H ,n( ds ) = [s] 1_2H ds. 

Let £t be the class of elementary functions on (0, T) x R d . For any if £ £® n , 

IMIW = d H [ f |/T-l /2 ([u] ff - 1/2 ^K, ■,..., U„,-)(0)(S)| 2 

JK" d J(0,T)" 

AH,n(rfs)/i(t^i).../i(t^ n ), (16) 

where| = (a,...,^)eM" d . 

The following theorem is the multidimensional analogue of Theorem 2.1 of 

0. 



Theorem 2.2 Let (0,T)™ 3t^ <p(ti, ■, . . . ,t n , •) £ S'(]R™ rf ) 6e a deterministic 
function such that !Ff(ti, •,•••, tn") *s a {unction on M. nd , for all t £ (0, T) n . 

Suppose that: 
(i) the function t h-> Ttp(t u ;...,t n , •)(£) belongs toU c {0,T)® n for all £ £ R nd ; 
(ii) the function (t, £) H ► F<p(ti, •,...,£„, ■)(£) is measurable on (0, T) n x K rad ; 

r«; x ( o,T)»nr=i iK>ajM K - i/2 [u-s] ff - 3 / 2 i^( U1 , •,..., Un ,-)(Didu < ^ 

/or aH (s,0 G (0,T) n x M» rf (or F<p(ai, -,...,*„, ■)(€) > /or aH (s,£)J. 
// 



It ■■= a 7 fj / J"y(ui, -,..., ■u„,-)(|)J"<^(wi, -,..., iw)(€) 

JK' ld J(0,T) 2 ™ 
n 

JJ |u< - ^| 2ir - 2 dudv M (dei) • • ■ M(d£n) < oo, (17) 

£/ien y G T-iP® n and \\ip\\ 2 ^ vSl - n — It- (By convention, we set <p(t\, •, ... ,i„, •) = 
if U > T for some i = 1, . . . , n.j 

Proof: The argument is similar to the one used in the proof of Theorem 2.1 of 
[2], being based on relations (JT3J) and ([15} above. We omit the details. D 



Remark 2.3 In our case, we apply Theorem 12.21 to the function 

(0,t) n 3t4 (f(ti, ■,..., t n ,-) = / n (<l, •,...,*„,♦,*, x). 

(We define f n (ti, •, ■ ■ ■ ,t n , •, t, x) to be if the relation t\ < ... < t n is not 
satisfied.) To see that this function satisfies hypothesis (i)-(iii) of Theorem 12.21 
we use ([9]) and the fact that: 

\TG(t,-)(0\<C T yt£[0,T}^£R d , (18) 

From here, we infer that the map t M> J-f n (ti, •,■••, t n , -,t, x) belongs to L^(0, T)® n , 
which is included in U c {0,T)® n . 

Therefore, to show that f n (-,t,x) £ HV® n , it suffices to prove that (fTT|) 
holds. This will be done in Section [3] 

2.3 Malliavin Calculus 

In this subsection, we introduce the basic elements of the Malliavin calculus 
with respect to the isonormal Gaussian process W (see [10] for more details). 

We first introduce the multiple Wiener integral with respect to W. Let 
g be the a- field generated by {W(ip); ip £ UP}. By Theorem 1.1.1 of [TO], 
L 2 (tt,g,P) = ®%L UT n , where UT n be the ra-th Wiener chaos of W. Hence, 
every F £ L 2 (il,Q, P) admits the following Wiener chaos expansion: 

oo 

F = E J « F ' ( 19 ) 

n=0 



where J n is the projection on T-LVrn for n > 1. By convention, JqF = E(F). 

The definition of the multiple Wiener integral /„ is similar to the white noise 
case (see Subsection 1.1.2 of [TO])- More precisely, I n is a linear and continuous 
operator from T-fp® 11 onto HV n - For any / £ / KV® r \ we write 

In{f)= \ f(t 1 ,xi,...,t n ,x n )W(dt 1 ,dxx)...W(dt n ,dx n ), 

J(R + xM d ) n 

even if f n is not a function in (t\,Xi), . . . , (t n , x n ). Note that I n (f) = In(f) and 

E(I n (f)I n (g)) = n\{f,~9)nv^, (20) 

where / denotes the symmetrization of /, i.e. 

f(ti, Xi, . . ■ ,t n ,X n ) — — 2_^ f{tp(l)t x p(l)t ■ ■ ■ itp(n)i x p(n))- 

' pes n 

(Here S n denotes the set of all permutations of {1, . . . , n}.) 

Any random variable F £ L 2 (tt,Q,P) admits the decomposition: 

F = J2ln(fn) (21) 

for some /„ £ WP® n symmetric, and 

oo oo 

E\Ff = 2 E\In(fn)\ 2 = £ U\ H/nll^n. 
n=0 n=0 

(By convention, /q = S(F) and Jo(a;) = a; for all x £ R.) 

We now introduce the derivative operator. Let S be the class of smooth 
random variables of the form 

F = f(W(<p 1 ),...,W(<p n )), (22) 

where / £ C£°(R n ), cp, e ?#>, n > 1, and C 6 °°(IR") is the class of bounded 
C°°-functions on W 1 , whose partial derivatives are bounded. The Malliavin 
derivative of F of the form ([22]) is an "HP-valued random variable given by: 

DF - E J^w^ • ■ • ' w ^)Wi. 

i—l 

We endow S with the norm \\F\g 1:2 := E\F\ 2 + E\\DF\\ 2 HV . The operator D 
can be extended to the space B 1 ' 2 , the completion of S with respect to || • || D i,a. 
The following result is the analogue of Proposition 1.2.7 of [TU] . 

Proposition 2.4 Let F be a random variable given by \21\) . If F £ O 1,2 , then 

D.F = Y,nI n -i (/„(-,•)), 

n>l 

where • denotes the missing (t, x) variable. 



Proof: It is enough to assume that F = I n (/ n ), for some symmetric elementary 
function /„ of the form /„ = J2™, ...,i„=i a Ji...»„ 1 A il x...xA,„, where m > n, 
Ax, . . . , A m are pairwise-disjoint bounded Borel sets in R + x K d , and a^...^ = 
if any two of the indices i\, . . . , i n are equal. Then 

In = J2 a *i-*™ E lA e(n)X-* A P(in) ( 23 ) 

l<il<...<i„<m peS({ii. ..i„}) 

In(fn) = n! J! a 4l ..,:„^(A 4l )...VK(A,J, (24) 

l<ii <...<i n <m 

where <S({ii, . . . , i„}) denotes the set of all permutations of {i\, . . . , i„} and 
IU(A) = W{l A ). Using dMl, the fact that D{FG) = {DF)G + F(DG), and 
DW(ip) = <p for any ip £ HV, we infer that 

D.F = n\ Y. °<i-*- E H (*) ( II ^( A '- 

l<ii<...<i„<m j=l \k^j 

Using (|23|) . we obtain: 

n-1 

/„-!(/„(•,•)) = ^ a h ... in Yl k W Wil%ft)) 

l<ii<...<i„<m pe^({ii...in}) i =1 

n 

l<ii<...<i n <m i=l pe5({ii...i n }\{^}) fe ^ 

n 

= (n-i)i 53 o*x-i-£u,wn w ^)- 

l<ii<...<i Ti <m j — 1 fc^j 

D 

The divergence operator S is defined as the adjoint of the operator D. The 
domain of 5, denoted by Dom 5, is the set of u £ £ 2 (f2; HP) such that 

\E(DF, u) HV \ < c(E\F\ 2 ) 1/2 , VF £ B 1 ' 2 , 

where c is a constant depending on u. If u £ Dom 5, then <5(it) is the element 
of L 2 ($7) characterized by the following duality relation: 

E(FS{u))=E(DF,u) HV , VF £ D 1 ' 2 . (25) 

If u £ Dom S, we will use the notation 

S(u) = [ f u(t,x)W(St,Sx), 

JO JM d 

even if u is not a function in (£, x), and we say that 8{u) is the Skorohod integral 
of u with respect to W. 



10 



The next result gives an important calculus rule, which plays a crucial role in 
the present article. This rule states, in particular, that the Skorohod integral of 
a multiple Wiener integral of order n coincides with a multiple Wiener integral 
of order n + 1, i.e. 



f n (ti,xi,. . . ,t n ,x n ,t,x)W(dti,dxi) . . .W(dt n ,dx n ) W(St,Sx) 

f n {t x ,xi, ...,t n , x n , t, x)W(dt 1 ,dx 1 ) . . . W(dt n dx n )W(dt, dx). 

: + xR d )"+ 1 

(26) 

Proposition 2.5 Assume thatu G L 2 (Tt;'HV) has the Wiener chaos expansion: 

«(•) = J2 W»('»*))> ( 27 ) 

?i>0 

where • denotes the missing (t,x) -variable, ■ denotes the missing n variables 
(ti,Xi), . . . , (£„, x n ), and f n is symmetric and lies in %'P® 71 fin the first n vari- 
ables). Then u G Dom 5 if and only if the series X)n>o In+i(fn) converges in 
L 2 (Q), where f n is the symmetrization of f n in all n + 1 variables. In this case, 



*(u) = $^»+l(/n) = X}W/»). 



ra>0 n>0 



Remark 2.6 (a) If u(t, x) is a function in (t, x), relation (|2T[) is interpreted as 
follows: for any (t,x) G R+ x R d , 

u(t,x) = J2 7 «(/nM, x)) mL 2 (n). (28) 

n>0 

(b) I£u(t,x) is a distribution in (i,x), relation (f27|) is interpreted as follows: for 

any e C£°(M+ x M d ), 

(u(«), <£) = £) /„((/„(•, •),</>)) inL 2 (n). (29) 

n>0 

(c) If u(t, x) is a function in t and a distribution in x, relation (127[) is interpreted 
as follows: for any t > 0, G C§°(R d ), 

(tt(t, *), 0) = 53 /„((/„(-, t, *), 0)) in L 2 (ft), 

n>0 

where * denotes the missing x- variable. 

Proof: Using the same argument as in the proof of Proposition 1.3.7 of [TO], it 
suffices to prove that for any G — I n (g) with g G 7-CP® n symmetric, we have: 

E{DG,u) H v = E{I n {f n ^)G). (30) 
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Without loss of generality, we may assume that g is a function in all variables. 
By Proposition 12.41 DG is a function given by 

D Siy G = nI n - 1 (g(;s,y)) Vs>0,y€R d . (31) 

We consider separately the following three cases. 

Case 1. u(t,x) is a function in (t,x). Using ([251) . (|3Tj) . the orthogonality of 
the Wiener chaos spaces, and (|20l) . we obtain 

E(DG,u) nr = a H E u{t,x){D s>v G)\t - s\ 2H ~ 2 f(x - y)dxdydtds 

J(R + xR d ) 2 

= na H / E(I n -i (f n -i(-,t,x))I n -i(g(-,s,y))) 

J(R + xR d ) 2 

\t - s\ 2H ~ 2 f{x - y)dxdydtds 
= n(n-l)la H (fn-i{-,t,x),g(-,s,y))))- HV mn-D 

J(R + xR d ) 2 



t - s\ 2H - 2 f(x - y)dxdydtd, 



.s 



= n\ (f n -i, g)-uv®" = nl (fn-i,g)ur® n = E (In(fn-i)I n (g)), 

where for the second-last equality, we used the symmetry of g. This proves (|30[) . 
Case 2. u(t,x) is a distribution in (t, x) (in S'(R d+1 )). In this case, we 
regularize f n as follows. Let tp £ Cfi°(M. d+1 ) be such that ip > 0, the support 
of ip is included in (0,1) x {x € R d ; \x\ < 1} and / Rd +i ip(t,x)dtdx — 1. Let 
ip e (t,x) = e~ d ~ 1 ip(t/e,x/e) and / n , e (-,») := Ve * /«(-,•), where * denotes the 
convolution with respect to the missing (t, ir)-variable, denoted by •. Note that 
/n,e (•)£>#) is a function in (£, a;) (see p. 245 of [H]). Let 

U e (t,x) = ^ i'nC/n.eO, *,*))• 

n>0 

We claim that u £ = tp s * u. To see this, note that for any <f> £ C^°(M. + x M. d ), 

(fofe *«)(•), = K.)^ e *^) = 53 /„«/„(•,. )^ e *^)) 

ra>0 

= ^/„((/„, £ (-,.)») = ^(/„(/^(-,-)),0) = (u £ (-)», 
n>0 n>0 

where we used (|2^|) for the second equality, and the stochastic Fubini's theorem 
for the second-last equality. 

Applying the result of Case 1 to u e , we get: 

E(DG,u s ) nv = E(I n (f n _ he )G). (32) 

Relation ([317)) follows by letting e — > 0. On the left-hand side of (|3"2"|) . we have: 

SUtie-ul&p = a Htd [ dr\T\ l - 2H f ^dO\^u £ (T,0~^u(T,0\ 2 

JR JR d 

= a H4 [ dr\T\-( 2H -V [ M (dO|^w(r,0| 2 |^ E (r,0-l| 2 ^0, 

JR JR d 
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as e — > 0, by the Dominated Convergence Theorem, where T denotes the Fourier 
transform in the (t, x)-variable and an.d is a constant depending on H, d and fx. 
On the right-hand side of (1521) . we have: 

E\I n (fn-l,e) ~ In{fn-l)\ = ft!||/n-l,e — fn-l\\-HV» n = 
~ n—1 p 

n\a n Hid dn...dT n - l dTT[\T i \ 1 - 2H \T\ l - 2H M(dei)...A*(de„-i)Mde) 

|J 7(n) /n-i(n,6.---,rn-i^n-i,r,0| 2 |^V' e (r,0-l| 2 ->0. ase->0, 

by the Dominated Convergence Theorem, where J r ^ n ' denotes the Fourier trans- 
form in all n variables (ti,a;i), . . . , (t n -i,x n —i), (t,x). 

Case 3. u(t, x) is a function in t and a distribution in x. The argument is 
similar to Case 2, based on a regularization of u in space. We omit the details. 
□ 

We now return to our framework. 

Definition 2.7 We say that u = {u(t,x);t > 0,x G M d } is a solution of ([1]) 

i/ u(0, x) = 1 /or a// £ G M d , and /or any t > 0, x G R d : 

(i) E\u(t,x)\ 2 < oo; 

(ii) u(t,x) is Ft-measurable, where F t = a{W s (A);0 < s <t,A G Bb(R d )}; 

(Hi) the process v^' x ' := G(t — •, x — *)u belongs to Dom 5 and 

u{t,x) = l + 5(v it ' x) ). 

Here, ■ denotes the missing s-variable, * denotes the missing y-variable, and 
G(t — s,x — *)u(s,*) denotes the multiplication of the distribution G(t — s,x— *) 
with the function u(s, *), for any s G (0,i). 

The following result concludes our preliminary discussion. 

Theorem 2.8 Suppose that for any t > 0, x G JR d ,n > 1. f n (-,t,x) G UV® n , 
f n (-,t,x) being the kernels introduced in Subsection \2. 1\ Then equation {IJ) has 
a solution if and only if for any t > 0, x G lR d , 

the series ^2 n>l I n (f n {-,t,x)) converges in L 2 (fl). (33) 

In this case, the solution is given by: w(0, x) — 1 for all x G M. d , and 

u(t,x) = l + J2 I Mn(-,t,x)), for allt>0,xeR d . (34) 

n>l 

Proof: Let v^ t,x ' be given by Definition ^. 71 We claim that v^ t,x ^ has the Wiener 
chaos expansion: 

„(*.-)(.) = ^J n (/ n+1 (. )# ,t 1 a:)), (35) 

n>0 

where • denotes the missing (s, y)-variable. 
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From (|35p . by Proposition 12.51 it will follow that v^' x ' G Dom 5 if and only 
if the series J2 n >o In+i(fn+i(-,t, x)) converges in L 2 (tt), and in this case, 

*(„('.»)) = J2 T n+l(U + l(;t,x)) = u(t,x) - 1. 
n>0 

It remains to prove (|35|) . If d < 2, then G(t, x) is a function in cc, and (|35l) 
is clear, since for any s > and y €E IR d , 

w^(a,i/) - G(t-« J s-|/)X) / n(/n(- > * > w)) = X)Jn(G r (t-a J a;-i/)/„(-,*,i/)) 

n>0 n>0 

= 5Z ^(/rt+lC", «, 2/,*, ^))- 

n>0 

Suppose now that d > 3. Then G(£, a;) is a distribution in x. Recalling the 
interpretation given to (|35p in Remark 12.61 (c). we show that for any s € (0, t) 
and^G Cg°(R d ) fixed, 

(«^( S ,*),0) = ^/ n ((/„ +1 (-, S ,*,t,x),0)) ini 2 (0). (36) 

n>0 

Since i/ i,x )(s, *) is the product between the distribution G(t — s, x — *) and 
the function u(s, *), the action of v <<t ' x \s, *) on <j> is given by: 

<y M) (s,*),0) = (G(*-s,a;-*),^u(s,*)) = [^u(s,*)*G(t-s,*)](a;) 
= ^2\<f>J n (s,*)*G(t-s,*)](x), 

n>0 

where J n (s,y) = I n (fn(',s,y)), and we used (JM]). 

To prove (|36p . it suffices to show that for any n > 0, 

[^(S, *)*G(t-S, *)](aO=In((/n+l (',*,*,*, »), <£))• (37) 

Let G e be a regularization of G in space, i.e. G £ (s, *) = ^ e * G(s, *), where 
V> £ = e" d V(x/e), V S G£°(M d ), V > 0, supp ip C {x e R d ]\x\ < 1} and 
J Rd r/>(x)dx = 1. Then ^G £ (s, *)(0 = FG(s, *)(0^e(0- 

Since G e (s, y) is a function in y, for any n > 0, 

[0J„, e (s, *) *G e (t-S,*)](x) =I n ((f n +l,e{-,S,*,t, x), (/))), (38) 

where J n , 6 (s, y) = I n {fn,e{-, s, y)), 

f n ,e(h,xi,...,t n ,x n ,s,y) = G s (s-t n ,y-x n ) . . ■ G e (t2-ti,x 2 -Xi)l{ tl< ... <tn<s y, 
fn+i,e(h,xi,. . . , £„, a:„, s, y, t, a;) = G e (t-s, x-y)f n>e (ti,xi, . . . , t n , x n , s, y)l{ s< t}- 



Relation (|37j) follows by taking the limit as e — > in (|38|) . This is justified 
below. 
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On the left hand side of (|55|) . we use the fact that for any j/gK d , 

E\J n ,s(s,y)- J n (s,y)\ 2 ^0, as e -> 0. (39) 

To see this, note that: 
E\Jn,e(s,y) - J n (s,y)\ 2 = E|7„ (/„ j£ (•, s, y)) -/„(/„(-, s,2/))| 2 = n!||(/„, e - /n)(-,s,y)||«-p<> 
= «!<#/ fl |*i - «i| aff - a / Hfn, e -fn)(ti, ;■.., t n ,-,s,y)($) 

F{fn,s ~ /n)(si, •, . . . , S„, ■, S, 2/)(£)M(d£l) • • • K d £n)dtds, 

where £ = (£1, . . . , £„) and t = (i 1: . . . , t n ). Here, /„(•, s, y) is the symmetriza- 
tion of /„(•, s, y) in the first n variables (ti,Xi),..., (t n , X n ), and T denotes the 
Fourier transform with respect to the missing variables x\,...,X n . 
Using ^, one can prove that: 



F(fn,e ~ fn)(tl,; • • • ,*n, ;S, y){£) = Ffn(h, ',..., in, ', S, S/)(C)[^e(^p(l)) 



where p is a permutation such that t p (u < ■ ■ ■ < t p t n ) ■ Relation (|39[) follows 
by the Dominated Convergence theorem, since T^j e {^) — » 0. The application of 
this theorem is justified since \Fip E ({;)\ < 1 and ||/„(-, S)2/)||fi£p®n < oo. 
On the right hand side of (|38[). we use the fact that: 

E\I n ((f n+1>e (-,s, *, t, x), 4>)) - I n ((f n+1 (-,s, *, t, x), 4>))\ 2 -+ 0, as e ->• 0. (40) 

To see this, note that 

£|£i((/n+l,s(-, «,*,*> »)><£)) - I n {{f n+1 (-,S,*,t,x),<p))\ 2 = n!||0 e (-, S, *, x)\\wpS 



./(Us") 2 " 7_i JR" d 



" I (0,*)*" ?=! 



Tg £ (si, •,..., s n , •, s, i, a;)(£)/i(d£i) ■ • ■ n(d£ n )dtds, (41) 

where £ = (&, . . . ,£„), t = (* X) ...,£„) and 

3 £ (-,s,i,a;) := ((/„+i, e - /„+i)(-,s, *,*,«), 0). 

Here, the action of 4> is on the missing y- variable (denoted by *), f n +i(-, s,*,t,x) 
is the symmetrization of f n +i{-,s, *,t,x) in the first n variables (ti,xi), . . . , (t n ,x n ), 
and the Fourier transform is taken with respect to the missing variables x\ , . . . , x n . 
Note that 



Fg e {ti, •,..., t n , ■, s, t, x){£) = / J"(/„ + i jE -/ n+ i)(ii, ■,..., tn, -, s, *, t, x)(£, £,)T4>(£,)d£, 
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where the first Fourier transform under the integral is taken with respect to the 
n + 1 missing variables x\,. . . ,x n ,y. Using ([9]), one can prove that: 

•7 r (/n+i, e -/n+i)(*i ! ;---,t n , ; $, *, t, x){£, = Tf n {ti, -,..., t n , •, a, x)(£)fc e , p (£, C), 

where 

M& = e-* x TG(t-s, ■)& + ... + £„+£) [MCp(l)) ■ ■ ■ MUi) + ■ ■ ■ + Zp(n) + 0-M 

and p is the permutation for which t p i\\ < . . . < t p i n y Hence 

Fg e {tl, -,..., tn,-,S,t,x)(€) = Tf n (tl, •,..., tn,-,S,x)($) / fe e ,p(£, 0^H0 d ^ 

JR d 

=: J7„(ti, ■,..., i„,-, S ,a;)(0^e,p(0- (42) 

Since Fi/j e (£) — > 0, by the Dominated Convergence Theorem, it follows that 
Ke, P {€) ~> as e ^ 0. (To justify this, we use (JTHJl and \Ti> e {£)\ < 1.) 

Relation (|40l) follows from (J?!]) and (l4"2"j) . again by the Dominated Con- 
vergence Theorem, whose application is justified by the fact that |if £jP (£)| < 
2C t J Rd |J"0(OK = : C t,<f>, and ||/„(-, s,x)||^,®„ < oo. D 

Remark 2.9 Let Uo(t, x) = 1 and u n (t,x) = 1 + J2k=i Ik(fk(',t, x)) for n > 1. 

Let Vn = G(t — ■, x — *)u n for any n > 0. Usin 
one can show that <5(i% ) = u„+i(t, x) — 1, i.e 



Let i4 = G(i — •, x — *)u n for any n > 0. Using the same argument as above, 



u n+ i(t,x) = l+ / G(i - s,x- g/)y n (s, j/)W(£s, £?/), Vn > 0. 

Jo Js. d 

In other words, {u T i}„>o plays the role of the Picard's iteration sequence used 
in the case H — 1/2. 

3 Existence of the Solution 

In this section, we examine condition (|33[) . in the particular case when / is the 
Riesz kernel of order a g (0,d), i.e. /(x) = c a4 \x\-^~ a \ (i(d£) = \t\~ a d£. 

Our main result shows that a > d — 2 is a sufficient for ([3"3"| , and hence a 
sufficient condition for the existence of a solution to ([1]) (by Theorem |2.8|l . 

Note that a > d — 2 is also the sufficient condition for the existence of a 
solution to ([1]), in the case when H = 1/2 and the solution is interpreted using 
a martingale measure stochastic integral (see Theorem 5.1 of [3]). 

Due to the orthogonality of the Wiener chaos spaces and (|2T)|) . condition (|33l) 
is equivalent to: 

S(t) := 1 + Y / ri\\fn(-,t,x)\\ 2 nv0n <co, (43) 

n>l 

where /«(-, i, x) is the symmetrization of /„(-, i, x) in the first n variables (t\, xi), 
. . . , (t n ,x n ). In this case, E\u(t,x)\ 2 = S(t). 
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We begin with the calculation of \\f n (-,t 7 x)\\'^ iv ^ n . At the same time, this 
calculation will show that f n (-,t,x) £ fiT® 11 (see Remark [273]) . (By abuse of 
notation, we use || • || w -p®« , even if we do not know yet that /„(-, t, x) e WP® n .) 

Note that 

\\fn(;t,x)\\ 2 HV ^=a n H f TTl^-^-l 2 "- 2 / Ffn(ti,;...,t n ,-,t,x)(Z) 



3 = 1 



J"/ n (si, •,..., s„, ;t,x)(g)(j,(d€i) . . . (j,(d£ n )dtds, 
where £ = (£1, . . .,£„), t = (*i, . .. ,t„) and s = (s x , . . . , s„). Let 

gl n \-,t,x) := n\f n (ti, -,..., *„,-,*, z). 

Hence, 

^/ n (i 1; •, . . . ,i„, -,i,x)(0 = ^.9t n) (-, t, x)(€), 
and 



/»(• ,*,aOllw>»» = ri^ a H I fl |*i-*j| 2H - 2 ^ (n) (t,s)dt*, 



3=1 

where 



^ (n) (M):= / ^ n) (-,t,s)(€)^ n) (-,t,a;)(0/*(dei). .•/*(««). (44) 

We let 

a„(t) := (n!) 2 ||/„(- , i,^)!!^^®^. 

With this notation, relation (j43|) becomes: 

S(t) = ^ -^a„(t) < oo. (45) 

We proceed to the evaluation of a n (t), which relies on the evaluation of 
?/>(™)(t,s). Using relation ©, one can prove that: 

Tg[ n \;t,x)(£) = e -*tti+--«»)- ^ ^G( Ul ,-)(ep(i))^G(«2,-)(^(i)+e P (2)) 

pes„ 



■■■FGiUn,- )(£ p(1) +-.-+Cp(n))l{t p(1) <...<t p( „ ) }, (46) 

where S n is the set of all permutations of {1, . . . , n} and Uj — t p (j+i) — t p (j) for 
1 < j < n, with t p (n+i) = t. 

In the argument below, since there is no risk of confusion, we omit writing 
the variable (t,x) of g\. (•,£, x). 
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It is known that, for any t > 0, G(t, •) is a distribution with rapid decrease 
in S'(M. d ), whose Fourier transform is given by: (see e.g. [15] ) 

^(V)(0 = ^fp, ^^ d - 

The following central result was announced in the introduction, and will 
allow us to estimate -</>(™)(t,s). 

Lemma 3.1 Assume that d — 2 < a < d. Then, 

Sm ^P |g " r,r a dt < C a , d t a - d+2 , for anyt>0,r 1 e K d . 
Proof: Using the change of variable £' = t£, we obtain, 



i = t a - d + 2 1 ■• ••;;;':■ -' i^ -///i-./i'. 

We claim that: 



sin 2 (l£'l),„ 

If 12 



/(a):= / ^yP|o-e|- a de<Ca, d , Vaelf 

./«<< IsT 

To see this, we change the variable a — £ into £, and we write 

Ka)-/ ^!^ierde+/ ^^ir^M*)^). 

J|5|<1 I? — a l J|f|>l Is — a l 

For ii(a), we use the fact that |2HL£| < i f or anv x > 0. Hence 

h(a) < [ \t\~ a dZ = c d [ \- a+d - l dX = c d -^—. 

J\£\<1 Jo d - a 

For 12(a), we use the fact that 

(see p. 81 of [13). In our case, t = 1. Hence 

h{a)<AJ —± |£|-^<4su P / — i jaiei-de. 

Finally, we observe that a > d — 2 is equivalent to J Rd 1+ | f ^ |£|~"rf£ < 00, 
which in turn is equivalent to 

sup / r^ — d£r Q ^<°° 

(see (5.5) of [3]). D 

Based on the previous lemma, we estimate ip( n \t, s). 
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Lemma 3.2 Iff is the Riesz kernel of order a > d—2, then for any t , s € [0,i]™, 

V3 (n) (t,s)<Q d [/3(t)/3(s)](«- d + 2 )/ 2 , 

where (3(t) = n™=i(*p(j+i) - t P (j))> ^( s ) = II"=i( s <xG+i) - s <x(:0)> and f/le P er " 
mutations p and a of {1, ... ,n} are chosen such that 

t p (i) < i p (2) < • • • < t p {n) and s CT (i) < s CT ( 2 ) < . . . < s CT(n) , (47) 

««£/l t p („+l) = S<x(„+1) = t- 

Proof: By the Cauchy-Schwartz inequality, 

^ (n) (t,s) < ^ (n) (t,t) i/2 v; (n) (s,s) i/2 . 

Let Uj = £p(j+i) — t p (j) for j = 1, . . . , n. Using (|4"4")) and (|4*5)) . we obtain: 

^(t,t) = / |^ n) (OI 2 M^i)-..M^n) 

^(ux, ■ )(e P (i))| 2 ■ • ■ \FG{u n , ■ )(£ p(1) + • . . + e P (n))| V(#i) • ■ • M<£0 

|^G(«i, • )(^)| 2 . . . IJ-GK, • )(£ + . . . + ^)|V«) • • • MO 
f sin 2 (ui |fl|) sin 2 ( M2 |^+^|) sin 2 ( Un |& + ...+&|) ,,,,_«,,, ,,, 

i»- kii 2 ' &+&i a ■■' i£+.--+e;i 2 Kl1 "" lu ^"" tn ' 

where we used the change of variable £'■ = £ p(j ) , j = 1, . . . , n. 
We now use the change of variable 

%=£+••• + & i = l,...,n. 

The inverse transformation is: ^ — ?7i> Cj = Vj ~ Vj-i>3 = 2, . . . , n. We get 



,„ / sin (m ??i / sin (u 2 M) , I -a 

v w (t,t) = / <% — : — ,2 — ??i / d?? 2 — : — ,2 — \m-m\ ■■■ 

sin 2 (u„|77„|) 

dr l n i 12 \Vn-Vn-l\ ■ 

'■ mr 

Using Lemma [3.11 iterativelv. we obtain ^^"^ (t, t) < C™ d (u\ . . . u n ) a ~ d+2 . 
The result follows. □ 

Proposition 3.3 If f is the Riesz kernel of order a > d—2, then for any t > 
and n > 1, 

5 «W < C(t) n *_ d+2 , (48) 



(n!) 
where C(t) = C aAH t 2H+a - d+2 . 
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Proof: Let h = (a — d + 2)/(2H). As in the proof of Proposition 3.6 of pQ, 
using Lemma 13.21 and inequality (16) of pQ, we obtain that: 



2H 

h. 



5»(*) < Cl dM (n\) 2H (f [(t -«„)... (*o-.s,)]" ( /s 

\"'0<Si<S2<. ..<S„<t 

" C °."- ff * U(n(l + ft) + l), 

^- ,^m +n(l+h)2H j; 

For the last inequality above, we used the fact that for a > 0, r(an + 1) = 
CVi(n!) a , where C n is a constant such that A~™ < C„ < A™ for some A > 1. □ 
The existence of the solution is immediate. 

Proposition 3.4 If f is the Riesz kernel of order a > d — 2, then condition 
\3S\) holds, and consequently, equation (QP has a solution. 

Proof: As we mentioned earlier, condition (f3"3")l is equivalent to (j4~5^1 . which in 
turn is satisfied, since by Proposition [ 



s v = E ^(*) ^ E JL < °°- 

n>0 ' n>0 V ' ; 

The second statement follows by Theorem 12.81 □ 

Since C(t) is an increasing function in t, the previous argument shows that: 

S(t) = V i-a n (t) < C T < oo, We[0,T], (49) 

* — ' n! 

n>0 

for any T > 0, i.e. sup/ t ^wq T i xR d -E|w(£,;r)| 2 < oo for all T > 0. 

Remark 3.5 In the case of the heat equation, it was shown in pQ that, if / is 
the Riesz kernel of order a > d — AH, then for any t > and n > 1, 

5 n (t) < <?(*)" 



( n !)-(d-a)/2' 

where C(t) = Ca,d,Hi 2H_(rf_a)/2 . In this case, 05) holds if a > d - 2. 

Remark 3.6 Using the same method as above, one can prove that J2n>o ~h a n(t) < 
C T < oo for all t £ [0,T], where a n (t) = (n!) 2 ||/ n (-,£,:c)ll«-p®™- 



20 



4 Moments of the Solution 

In this section, we show that the solution is L 2 (Jl)-continuous and has uniformly 
bounded moments of order p > 1. With the obvious modifications, the results 
presented in this section remain valid for the heat equation (see Remark 13. 5[) . 

Let u(t, x) — J2 n >o Jn{t, x), where J n (t, x) is the projection of u(t, x) on the 
Wiener chaos T-UP n - By the orthogonality of the J„(i,x)'s, we have: 

n>0 

Note that 

E\J n (t,x)\ 2 =E\I n {f n (- ,t,x))\ 2 = -n n {t). (51) 

n\ 

It is known that, for any 1 < p < q < oo, the norms || • || p and || ■ \\ q 
are equivalent on any Wiener chaos 'HfPn, where || • || p denotes the norm in 
L p (fl). This is a consequence of the hypercontractivity property of the Ornstein- 
Uhlcnbeck semigroup (T t )t>o, defined by: 



T t F = J2 e ~ Ut J " F > FeL 2 (n), 



n>0 

where we denote by J n F the projection of F on the n-th Wiener chaos HVn- 
The property says that for any p > 1 and t > 0, 

\\T t F\\ q(t) < \\F\\ P , 

where q(i) — e 2t (p — 1) + 1 (see Theorem 1.4.1 of [ID])- Hence, for any 1 < p < 
q < oo and for any F G HV n , 

e- nt \\F\\ g = \\T t F\\ g < \\F\\ P , 

where t > is chosen such that q = e 2t (p — 1) + 1. In particular, for any p > 2 
and for any F £ WFn, 

||^|| P < e "lF|| 2 = (p-l)»/ 2 ||F|| 2 , (52) 

where t > is chosen such that p — e 2t + 1 . 

Applying these results in our case, we obtain the following result. 

Theorem 4.1 Let f be the Riesz kernel of order a > d — 2 and u be the solution 
of (QP. Then u is L 2 (fl)- continuous, and for any p > 1, T > 

sup sup E\u(t, x\ p < oo. (53) 

t<T X £R d 

Proof: We apply J52]) for F = J n (t,x) G HV n . Using (J5TJ) , we obtain that: 

/ -, \ 1/2 

\\Jn(t,x)\\ p < (p~l)^ 2 \\J n (t,x)\\ 2 = {p- I)"' 2 \-fn{t)J ■ 
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Using (|4"5)l , we obtain that: 

£ \\Ut,x)\\ p < J> - 1)" /2 WTT^m V /2 < oo. 

n>0 n>0 ^ V '' J 

Since a(t) does not depend on x and C(i) is an increasing function of t, we 
have: for any T > 0, 

y^ sup sup || J n (t,x)\\ p < C T .p < oo- (54) 

n > t<T s6 «- 

From here, we conclude that for any (t, x) £ R+ xM. d , the sequence {u n (t, x) — 
Y^k=o Jk(t, x), n > 0} is Cauchy in L P (Q), since 

n 

\\u n (t,x) — u m (t, x)\\ p < \, ll>^fc(i, a;)||p — > 0, as n, to — > oo,n > to. 

fc= 771+1 

Therefore, there exists a random variable u(t, x) £ i p ($7) such that u n (t,x) — > 
v(t,x) in L p (fl). But u n (t,x) — > u(t,x) in F 2 ($7), and hence u(t,x) = v(t,x) 
a.s. Using (|54l) . one can show that 



" 



( (t, a?) -> w(i, x) in L p (fi), uniformly in (£, x) G [0, T] x R d (55) 



and ||w„(i,x)|| p < Y2=o \\ J k(t,x)\\ p < C T , P for all (t,x) £ [0,T] x R d ,n > 0. 
Taking n — > oo, we obtain (|53[) . 

By LemmaS2]below, J n is F 2 (f2)-continuous. Hence w n is F 2 (f2)-continuous. 
Due to ([55]) . it follows that u is L 2 (J7)-continuous. □ 

Lemma 4.2 aj For any n > 1 and t > 0, 

E\J n (t + h,x) — J n (t, x)\ — > as h — > 0, uniformly in x £ M. d . 

b) For any n > 1, t > 0, x £ R d , 

E\J n (t,x + z) - J n (t,x)\ 2 ->-0 as \z\ ->0,ze M d . 

Proof: a) Suppose that h £ [0, 1]. (The case h < is similar.) Then, 

£|J„(i + /i,a:)- >/n(i,z)| 2 = E\I n (f n (-,t + h,x)-f n (-,t,x))\ 2 

= »! ||/n(-,* + ^,a:)-/n(-,*)a;)||^73®» 

n! 



where 



£i(t,ft) := (nl) 2 \\f n (',t + h,x)l m » -/„(-, t,x)f nv ^ (56) 

E 2 {t,h) := (n!) 2 ||/„(-,i + /i,x)l [O!f+/i] „\ [ot] „|| 2 ^ 7 r )0 „. (57) 
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We treat Ei(t, h) first. Note that 

« n 

E x {t, h)=a n H IT \tj - Sj \ 2H - 2 ^\t, s)dtds, (58) 

•Ao,*] 2 " y=i 



where 



V>£ n) (M) = / J-(9 t (n H^,x+/ l )- fft ( " ) (^,x))(OJ-( ff |™ ) (^^ + / l )- 5 i n) (-,t,x))(OM(rf6).--M(^) 



By the Cauchy-Schwartz inequality, ^(M) < ^ n) (t, t) 1 / 2 • ^(s^) 1 / 2 . 

To evaluate ip^ (t,t), we use ()46[) . denoting Uj = t p (j+i) — ^p(j)^ when 
< *p(i) < • • • < t p („) < i p (n+i) = i: 



|^G(«i, • )(^ p( i))| 2 ■ • • I-FGK-!, . )(^ (1) + . . . + £ p( „_ 1} )| 2 
\T[G(u n + h,-)- G{u ni • )](e p( i) + . . . + e p (n))| 2 |«(^i) • • • A*(#»). 
Proceeding as in the evaluation of i/>*( n )(t, t), we obtain, 
,(«)/+ .n /" , sin 2 (Mi|??i|) y sin 2 (u 2 |7? 2 |) | ,_ Q 

[ sin(K + fr)Kl) - sinKKi)] 2 _ _ a _ 

By the Dominated Convergence Theorem, 

^ n) (t,t)^0 as/i->0. 
The application of this theorem is justified, since 

\ S m((t + h)\S\)- S m(t\l;\)\ < ( 4 V /2 



iei -vi + iei 

for all £ e R d ,i > 0,h G [0,1] (see p.4 of Erratum of [4). The fact that 
E\(t, h) — > follows by applying the Dominated Convergence Theorem in (|58p. 
We now treat E 2 (i, M- Let ^ = [0. * + h] n \[0, t] n . We have 

„ n 

E 2 {t,h) = a n H T\\t - S:l \ 2H - 2 l A (t)l A (s) 7 i n \t,s)dtds, (60) 

where 



7 i n) (t,s)= / Tg { t n \- ,t + h,x)(£)Tgi n) (- ,t + /i^)(^6)-M4)- 
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By the Cauchy-Schwartz inequality, 7£ n) (t,s) < 7^ l) (t,t) 1 / 2 7 < ( n) (s,s) 1 /2. To 
evaluate 7^ (t,t), we use again (|46|): 



» 



7,r(t,t) 



\Fgi n \-,t + h,x)(t)\ 2 »(<%!)...»((%„) 



sin 2 (Mi|??i|) 






sin 2 (u 2 |'72| 






l??2 -??i| 



(61) 



where ttj = ^(j+i) — ^o(j')- By Lemma 13. 11 



» 



7r ; (t,t)<CS )<J [u 1 ...« n _i(« n + /i)] 



a-d+2 



and hence 



» 



7 r(t,s)<c, d 



(u n + h)(v n + h) YJ_ 



3=1 



(a-d+2)/2 



where Vj = s a{j+1) - s a{j) and s ct(1) < . . . < s ff(n) < s ff( „ +1) = t. 
By inequality (16) in [T], 



E 2 (t,h)<b 2 H n Cl, 



( 



[0,t+h] 



U(t) 



ra-1 



II( t PO'+ 1 ) ~ f P0))( f + h - *P(n)) 



2ff 



r/t 



where 5 = (a - d + 2)/(2H), Using the fact that 

A = IJ {(*!)•••>*"); °< */>(l) < ••• <*p(n-l) <*n aildt n £ [t,t + /l]}, 



we obtain that: 

£ 2 M) < &£<%„ 
where 



-i 2H 



t+h 



n-l 

JJ («,-+! -t i ) ff dti...dt„_i 

0<ti<...<t„_i<t n J=1 

r(i + <s)™ 



r((n-l)(l + <J) + l) n 



x(n-l)(l+<5) 
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(see Lemma 3.5 of pQ). We obtain 
E 2 {t,h) < b%C2 d T(l+5f Hn 



-l2H 

t+h 



fa 12! f (t + h-t ^t^-W+^dt 

r(( n -l)(l + 6) + l)J t [+ n) n n 



t+h 



= cz AH {t + i)mn-m + s) _ _ L_ h 2Hv + i) _> , as /wo. 

b) We have: 



S|J n (i,x)-J„(i,2/)| 2 = £?!/„(/„(■, i, x)) -/„(/„(-, t,y))| 2 = - f ^ 3 (t,a:,l/) J 



1 
where 



£3(^,2/) := (ri!) 2 ||/„(-, t, a) -/»(-,*, J/)||^» n (62) 

W n^-^| 2H - 2 ^(t,s)dtds (63) 



and 



^™>(M) = / ^(flr(-*»*)-5r(-.*.i/))(€)J 7 0;i" , (-t J »)-fli B) (- ,t,i/))(^(€)-^n). 



By the Cauchy-Schwartz inequality ip^(t,s) < ipi™y(t, t) 1 / 2 ^i"^(s,s) 1 / 2 . 
To evaluate ^4™y(t,t), we use (|15)) : 

^(t,t) - / |^(ff t (n) (- ) t,x)- fft (n) (-,i, 2 /))(0|V^l)---M(^n) 

sin 2 (Mi|^i|) /" sin 2 (w 2 |77 2 |) 
^ T~i2 W / d m 7—T2 »72-»7i ••• 

I ml .V M 2 

d ^ sin 2 K|77 n |) |7?w _ ??n _ ira|1 _ e -»,„.fo-)|2. (64) 

By the Dominated Convergence Theorem, E 3 (t, x, y) — > as |a; — y\ — ► 0. D 

5 Holder Continuity 

In this section, we obtain some bounds for the p-th moments of the solution, 
from which we infer that the solution has a 7-Holder continuous modification, 
with0< 7 < 2=|±2 # 

If / is the Riesz kernel of order a > d — 2, then for any 1 > /3 > (d — a)/2, 
/ltd n + ifpV < °°' which is equivalent to 

SU P / n l^-lira < 00. (65) 
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(see (7.26) of 0). 

By Proposition 7.4 of [3], the fundamental solution G of the wave equation 
satisfies hypothesis (H3)-(H5) of [3], for any < 7* < 1 -(3 < 2=|±2, j = 1, 2, 3. 
This fact is used in the proof of the next result. 

Theorem 5.1 Let f be the Riesz kernel of order a > d — 1 and u be the solution 
of $7$). Then for any p > 2, T > and K C K d compact, 

E\u{t + h,x)~u{t,x)\ p < C\h\ phl ^ 2+H)] , VtE[0,T},Vx£m. d ,Vh£R,t + h£[0,T}, 
E\u(t,x + z)-u(t,x)\ p < C|z| P73 , \ft>0,Vx€K,Vz€R d ,x + z€K 



for any < 7j < " 2 + 2 ; where C is a constant which depends on a,d, H,p,T. 
In particular, {u(t, x); (t, x) € [0,T] x K} has a modification which is a.s. 
jointly j-Holder continuous in time and space, for any 7 £ (0, a ~^ +2 ). 

Proof: We first treat the time increments. By Minkowski's inequality and (|52"j) . 



\\u(t + h.x) - u(t,x)\\p = || 22^ n ^ t + ^' x ) _ Jn{t,x))\\ p 

n>0 

< J2 W J n(t + h,x)- J n {t, X)\\ p <J2(P~ i)" 72 !! J n(t + h,x)- J n (t, x)\\ 2 
n>0 n>0 

= Y,(p-l) n/2 {^(Ei(t,h)+E 2 (t,h))X , (66) 

n>0 *• ' ' 

where Ei(t, h) and E 2 {t, h) are given by ([56]) . respectively (|57|) . 

To estimate -Ei(t, h), we use (|58|) and (|59|) . The inner integral in (|59|) is 

|J-GK + ft, • )(£ + r/^O - JG( U „, • )(£ + ?7n-i)| 2 M(^)- 

This integral is bounded by Ch 2 " 11 for some < 71 < °^+ 2 , due to (H3). The 
remaining (n— l)-fold integral in (|59|) is bounded aboveby C™"^ 1 (mi . . . u n ^i) a ~ d+2 , 
by Lemma |3. 11 Hence, 

</4 n) (t,s) < Ch 2 ^ci; d \ Ul . ..u n ^ a - d+2 ^ 2 ( Vl . ..v n -i)^- d+2 ^ 2 , 

where Uj = t p (j+i) — t p (j) and Vj = So-q+i) — s a(j)- By inequality (16) of pQ, 

n-1 



2 if 



£a(*,ft) < Ch^C^b^inl f f Y[(t j -t j -i) s dti...dt n -idt n 

\ JO J0<ti<...<t n -i<t„ j =1 

/ ! /"* \ 2 " 

- C/i 27l C™" 1 6^rd + (5 , ) 2H " I — / i (n ~ 1)(1+,5) d£ I 

- ^ ^m^UIH-0) ^r((n-l)(l + 5) + l)y *» * r V ' 

Tin(2H+a-d+2) 

< Cft J1 C™ AH ,N Q _ d+2 , (67) 
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where 6= (a - d + 2)/(2H). 

To estimate E2(t, h), we use (|60|) and (|6"Tj) . The inner integral in (|6Tj) is 



|^GK + /»,-)« + ry„_i)| 2 /x(dO < C(«„ + /i) 2 ^, 
for some < 72 < a ~^+ 2 ; by (H4). Using the same ideas as above, we get: 

Mt,h) < cc ldtH ( j n m W - f\t + h - t n) t^^ dt 



< Ch«-»+*C2 tdtH 2 . (68) 



r((n-i)(i + *) + i)y t 

rpn(2H+a-d+2) 

n!) 
From ([66 ]l . (|6? j) and ([68 ]) . we get: 

(p - l)™/2(7"/2 T n(2H+a-d+2)/2 



\\u(t+h,x)-u(t,x)\\ p < Ofc^^+^X; fn n(«-, + 3)/ 2 

n>0 ^ '-' 

=: h^ h ^ +H ^C{a,d,H,p,T). 

We now treat the spatial increments. As above, we obtain: 

/ 1 x 1/2 

||u(t,x + *)-u(t )a 0|| p < ^(p-1)™/ 2 (-E 3 (t,x,x + z) 



x n! 

n>0 x 



where E 3 (t, x, y) is given by (|52"j) . To estimate E^x, x + z) we use (|6"5t and (|M)l . 
Using the fact that J-G(u, ■ — z)(£) — e~ l ^' z FG(u, •)(£), we see that the inner 
integral in (IM1) is: 

I JG( Un , • - z)(£ + Vn-i) - FG{u n , ■ )(£ + Vn-i)\ 2 KdO < C\z\ 2 ^, 

for some < 73 < °~^+ 2 , by (H5). The rest of the proof is the same as above. 
The final statement follows by a version of Kolmogorov's criterion for multi- 
parameter processes (see e.g. Problem 2.9 of [8]). □ 

6 Malliavin differentiability of the solution 

In this section, we show that u(t,x) is Malliavin differentiable of any order. 
When d < 2, we show that the Malliavin derivative of the solution satisfies a 
certain integral equation. These results are valid for the heat equation in any 
dimension d. 

Recall that if F is a smooth random variable of the form (|22|) . the iterated 
Malliavin derivative D k F is an 'H'P® -valued random variable, defined by: 

° kF = E a {W{vi),...,W{<p n ))<p ix ®...®ip ik . 

OXi-. . . . Xj, 



41,.. l k = 1 
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The space lD> fc ' p is the completion of the space S of smooth random variables, 
with respect to the the norm || • ||d*.p defined by: 

k 
\\Ff BkiP = E\FF +J2\\D j F\f Hr ^. (69) 

Using a known criterion (see e.g. p. 28 of [10), it follows that u(t, x) £ O fe ' 2 , 
since by ([ST]) and (|48|) . we have: 

y: n k E\ut, x) \^y: - k ^n(t) < e |S§& < ~. 

n>l ra>l ' n>l V '' 

Next, we show that u(t,x) € D fe - P for all k > l,p > 1. 

Let F g L p (f2) be such that £> fe F exists. By Meyer's inequalities (Theorem 
1.5.1 of [rU]), for any p > 1, 

FIIZ^FII^,^ < oo if and only if E|C fe F| p < oo, (70) 

where C fc : Dom C k C L 2 (fi) -» L 2 (f2) is the operator defined by 

n>l 

and Dom C k = {F e L 2 (fl); J2 n >i n k E\J n F\ 2 < oo} = D fc < 2 for any fc > 1. 
By Minkowski's inequality and (|52p . we have: 

r fe F|| p <En fe / 2 ||J„F|| p <En fc / 2 (p-l)"/ 2 ||J„F|| 2 . (71) 

n>l n>l 

Combining ([55]), ([701) and ^7TJ), we infer that ||.F|| D *,* < oo (i.e. F € © fc ' p ), 
if 

En fe / 2 (p-l)"/ 2 ||J n F|| 2 <oo. (72) 

n>l 

Applying this in our case, we obtain the following result. 

Proposition 6.1 Let f be the Riesz kernel of order a > d — 2 and u be the 

solution of (pp. Then u(t,x) £ H) k,p for all k > 1 and p > 1. 



Proof: We verify ([72]). By (|5T) and (|48[). we have: 

E^(l- l) B/2 ||Jn(t,x)|| 2 = ^ n l/2 ^l)" /2 (ia,(t) N 

n>l n>l ^ 

[2 fe (p-l)C(£)]"/ 2 

- Z^ f„lMa-d+3)/2 ^ °°- 

n>l v ' y 

□ 
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In the final part of this section, we show that the Malliavin derivative 
Du(t,x) satisfies a certain integral equation. For this, we assume that G(t,x) 
is a function in x (i.e. d < 2). 

Recall that u satisfies the integral equation ([3]). Intuitively, using the com- 
mutativity between the operators D and S, the derivative Du should satisfy: 

Du(t,x) = G(t--,x-*)u+ I [ G(t-s,x-y)Du{s,y)W(6s,Sy), (73) 

JO JR d 

where • denotes the missing r variable and * denotes the missing z variable. 

The integrand of the stochastic integral above is an HV (8)WP -valued random 
variable, and the integral needs to be defined as an 'H'P-valued random variable. 

For this reason, we introduce a Hilbert-space- valued Skorohod integral. 

If A is an arbitrary Hilbert space, we let S(A) be the class of smooth A- 
valued random variables F — X)i=i Fj v j> with Fj £ S,Vj £ A,m > 1. The 
Malliavin derivative of such F is defined as DF = J2T=i DF <g> Vj. We denote 
by D 1,2 (^l) the completion of S(A) with respect to the norm || • Hb 12 ^), where 



pi.2^) := E\\F\\ A + E\\DF\\ HV9A . 



Similarly to the case A = M (considered in Subsection 12. 3p . we let S* be the 
adjoint of the operator D. The domain of 5* , denoted by Dom S* , is the set of 
U e L 2 (fl; KP ® A) such that: 

\E(DF, U)uv®a\ < c(E\\ F\\ A ) x /\ VF G B^ 2 (A), 

where c is a constant depending on U. If U £ Dom S* , then S*(U) is the element 
of L 2 (£l;A) characterized by the following duality relation: 

E{F,5*(U)) a = E{DF,U)hv®a, VFeB 1 ' 2 ^). (74) 

If U £ Dom 5*, we use the notation 

S*(U)= I I U(t,x)W{S*t,S*x) 7 

JO JR d 

even if U(t, x) is not a function in (t, x), and we say that S*(U) is the A- valued 
Skorohod integral of U with respect to W . 

Similar to the case A = M, we have the following result. 

Proposition 6.2 Let U, V £ B^ 2 (nV ® A). Then 

(OO 
J2 D e * (U, e 3 a k ) HV9A Di (V, e % ® a k ) U v®A 

where {ei)i>\, (a.k)k>i are complete orthonormal systems in T-LV ', respectively 
A. Consequently, ifUe O l ' 2 {UV ® A), then U £ Dom S* and 

-^11^(^)11^ < E W\\nv®A + e \\ d U\\hv®uv®a = \\ u \\w^{-hv®a)- ( 75 ) 
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Proof: The proof is similar to Proposition 1.3.1 of [TU]. For this, one needs to 
revisit the basic rules of the Malliavin calculus. We omit the details, but we list 
below these rules, which apply to any isonormal Gaussian process {W (h)} hen'- 

1) For any F E S{A),hE H,V E A, 

E(DF, h <g> v) n ® A = E((F, v) A W{h)). 

2) For any F E S(A), G E S, h £H, V E A, 

E(G(DF, h ® v) H9A ) = -E((F, v) A (DG, h) H ) + E((F, v) A GW{h)). 

3) If U = YZLi Fj(hj ® Vj ) E S(H <g> A) for some F j E S, h 3 E H, Vj E A., 
then U E Dom 6* and 

m m 

6*(U) = J2 FjW(hj) - ^(DF^h^nVj. 

3 = 1 3=1 

4) For any U = ~£%i Fj{h 3 ® Vj) £ S{H ® A) and h E H, v E A, 

D h ^(S*(U)) ={U,h® v) n ®A + (6*(D h u),v) A , 

wheTeD h ® v (6*(U)) = (D(6*(U)),h<S>v) n ® A and D h U = Y™Li(D h Fj)(hj®Vj)- 
Here D h F = (DF,h) H . D 

In what follows, we let 6* be the operator corresponding to the case A = HV. 
We begin with some preliminary results. 

Lemma 6.3 For any t > 0,x E R d , Du n (t,x) -> Du(t,x) in I? '(ty'H'P) , and 

C ( t ] := sup E\\Du(t, x^-p <oo, forallT>0. (76) 

(t,x)e[0,T]xR d 

Proof: Using Lemma 1.2.3 of [10], it suffices to prove that: 
sup sup E\\Du n (t, x)\\-wp < co. 

">1 (t,x)e[0,T]xM d 

By Proposition H2D D rtZ u n (t,x) = YX=i kI k-i{fk{- ,r,z,t,x)). Using the 
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orthogonality of the Wiener chaos spaces, ([20)) and (|4l))) , we get: for any t £ [0,T] 

E\\Du n (t,x)\\ 2 nv = a H [ f \r-r'\ 2H - 2 f{z-z')(y^kI k ^{f k {- ,r,z,t,x)) 

J(o,ty Jr™ \^[ , 

Ij2 ll l-M- ,r',z',t,x))\ dzdz'drdr' 

= a H J2k 2 {k-l)\f f \ r -r'\ 2H - 2 f(z-z') 
%T X J(Q,ty Js** 

(fk(- ,r,z,t,x),fk(- ,r' \z ',t,x)) nv ®(k-i)dzdz'drdr' 

n n 

— / J kk\ \\fk{' )t)%)\\wp®ii f? / J kkl ||/fe(- ,t,x)\\ HV(Sk 

fc=l fc=l 

n _. n 1 

fc=i v '' fc=i v '" 

where fk(-,r,z,t,x) denotes the symmetrization of fk(-,r,z,t,x) with respect 
to the first k — 1 variables. For the first inequality above, we used the fact that 
the H/H^-pgm < WfWuv®"- f° r an y / G WV® n . For the last inequality, we used 
Remark EH □ 

Remark 6.4 Using Proposition ^. 41 iterativelv. we obtain that: 

D fr,w),(r,z) U ( t ^ X ) = ^2 n ( n ~ 1 ) / "-2(/n(', T, W, r, Z, t, x)), 
n>2 

for any (r, w) £ (0, t) x R d and (r, z) £ (0, t) x R d . Hence, similarly to (1751) . one 
can show that: 

C^ 2) := sup E\\D 2 u{t,x)\\ 2 nv ^ nv < oo, for all T > 0. (77) 

(t,x)e[0,T]xR d 

Lemma 6.5 For any t > 0, x £ W 1 , the process U^' x ' defined by 

jj(t.x) ={J7(*.») («,|,) :=G(t-s,.x- 2 /)L»u(s,y);s>0,yeK d } (78) 

belongs to Dom (5*. 

Proof: By Proposition [6~2l it suffices to show that U^ £ n l > 2 (WP ® WP), 
i.e. ||£^'' x )||B 1 . 2 (H'Pig>w) < °°. Note that 

iirr(t,a;)||2 _ tp\\ rr(t,x) ||2 , p|| r)TT(t,x) ||2 
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By the Cauchy-Schwartz inequality, (175)) and ([77)1 , 



E\\U\\ 2 nvm = a H [ [ G(t-s,x-y)G(t-s',x-y)E(Du(s,y),Du(s',y% 

Jo Jo 



\s — s 



J\2H-2 



f{y-y')dydy'dsds' 



< Cl x '\\G(t-;x-*W nv <cx>, 

E\\DU\\ 2 nv ® 3 = an f [ G{t - s,x - y)G{t- s' ,x - y)E(D 2 u(s,y),D 2 u(s',y')) npcS 2 
Jo Jo 

\s-s'\ 2H - 2 f(y-y')dydy'd S d S ' 

< C( 2) \\G(t--,x-*)\\ 2 nv <oo. 

D 

Using the same argument as above, one can show that the process C/„ 
defined by C/„ ,x (s, y) := G(t — s, x — y)Du n (s,y), belongs to Dom 6*. 

The next result shows that the sequence {Du n (t, x)} n >o satisfies a recurrence 
relation. 

Proposition 6.6 For any t > Q,x € K d and n > 1, 

DUn(t,x) = G(t-;X-*)fln-l+ [ / u£$(a, y) W(8* s, 6*y) 

Jo Jr- 1 

mL 2 {n;HV). 

Proof: Step 1. By induction on n, we show that for any r e (0,i), z G M d , 
D r ^ z u n (t,x) = Y,k=i Ak ' wnere A k = Xh=i ^fe i 

-4fe = / G(t - Sk-i, x- yk-i) ...G(si - r,yi - z)W(ds 1 ,dy 1 ) . . .W(ds k -i,dy k -i) 

Jr<si 

A k l) = l G(t-Sk-i,x-yk-i).-.G(si-r,yi-z)G(r-Si-i,z-yi-i) 

Jsi—i <r<Si 

■ ■■G(s 2 - si,y 2 -yi)W(ds]_,dyi) . . . W(dsk-i,dyk-i), i = 2,...,k-l 
A k > = / G(t-r,x- z)G(r- s k -i,z-yk-i)---G(s2- si,y 2 -yi) 

J s k -i<r 

W(ds 1 ,dy 1 ) . . . W(ds k -i,dy k -i), 

and the integrals above are taken over the set {0 < s% < . . . < Sk-i < t} x M. nd . 
By definition, u n (t,x) = u n -i(t,x) + I„(f n (-,t,x)). Using Proposition 12.41 
and the induction hypothesis, we obtain that: 



D rtZ u„(t,x) = ^A k +nl n -i(f n (-,r,z,t,x)). 
fc=i 

Note that n/„_i (/„(-, r, z, t, x)) — A n , the n integrals A" , . . . , An correspond- 
ing to the n possible locations of r, compared with the variables Si < . . . < S n _i. 
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Step 2. We prove that for every r e (0, i), z e R d 
D r>z u„(t,x) = G(t-r,x-z)u n -i(r,z)+ / / G(t-s,x-y)D r>z u n -i(s,y)W(5s,5y). 

Jo JR d 

We use the expression of D rz u n (t, x) obtained in Step 1. Note that the terms 

A\ ,A 2 , • • ■ , An have the common factor G(t — r,x — z). A quick calculation 
shows that the sum of these terms is G(t — r,x — z)u„_i(r, z). 

For the remaining terms, we change the names of the variables of integration, 
so that G(t — s,x — y) becomes a common factor. More precisely, we call (s, y) 
the variable (sk-i, 2/fc-i) in A^ , for any i = 1, . . . , k — 1 and k — 2, . . . , n. The 
sum of these terms turns out to be 

/ G(t- s,x - y)D r . z u n - 1 (s,y)W(ds,dy), 

JR d 

using the expression of Z?,- iZ u„_i(s, y) obtained in Step 1. Finally, by (|2"5)l . we 
can replace the integral W(ds, dy) by an integral W(Ss, Sy). 

Step 3. We show that the process 6*(U%fJ) coincides (in L 2 (n,HV)) with 
the process V^l* , defined by 

V£.?(r,z) ■= [ [ G(t-s,x-y)D r!Z u n ^(s,y)W(6s,6y) 

Jo J«. d 

= 6(G(t- -,x- *)D r>z U n -l). 
By the duality relation (jTl)) . it suffices to prove that 

E(F, v£$) HV = E(DF, Ujffhnrwp, VF e U ll2 («P). (79) 

Without loss of generality, we may assume that F is smooth, i.e. F = F^ip 
with F € S and </? e UV . Then DF = £>F ® ip and 

£(DF, UJffiHrvHV =a H [ E(DF , G(t -*,x- -)D riZ u n ) nv 

J(M + xM d ) 2 

(p(r',z')\r - r'\ 2H - 2 f{z - z')dzdz'drdr' 
= a H J E(FQV^{r,z))tp(r',z')\r - r'\ 2H - 2 f(z - z')dzdz' drdr' 

J(R + xS. d ) 2 

= E{F,V { - t ' x) )nv- 

Note that for the second-last equality above, we used the duality relation (|25p 
(for the operator S), whereas for the first and last equality we used Fubini's 
theorem. This shows (f79|). and concludes the proof. □ 

The following result gives the precise meaning of relation (I75|) . 
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Theorem 6.7 Let f be the Riesz kernel of order a > d— 2 and u be a solution 
of (Qp. For any t>0,X<= R d , let U (t ^ be defined by |73p. Then, 

Du(t,x)=G{t--,x-*)u+ [ [ U {t ' x) {s,y)W(S*s,S*y) 

inL 2 (n;HV). 

Proof: By the duality relation (|74p . it suffices to prove that, 

E(Du(t, x)-G(t--,x- *)«, F) uv = E(DF, U^) nv ® nv , (80) 
for any F G B^ 2 (HV). By Proposition HH for any F G B^ 2 (nV), 

E(Du n (t,x) - G(t --,x- *)u n - X ,F) nv = E{DF,U { ^)uv®uv- (81) 



Relation (|80|) is obtained by taking n — > oo in (J8TJ) . We justify this below. 
On the right-hand side of (j8Tj) , we use the duality relation ([74]) . the Cauchy- 
Schwartz inequality, and ([75]) : 

E(DF, Utl ~ U^) Hvmv - E(S*(ut1 - U^ x) ), F) uv 

< (E\\5*(U^ U^)\ 2 nv) l,2 {E\\F\\ 2 nv) l/2 

< ll^n-i -^ \\v>i,2Chv®hv)( e \\ f \\hv) ■ 

To show that ||17„_li — U ( - t ' x ' ) \\ 2 llli 2(- H 'P®'H.v) — > as n — » oo, we use the same 
argument as in Lemma 16.51 

For the first term on the left-hand side of (f8"Tj) . by the Cauchy-Schwartz 
inequality and Lemma 16.31 



E{Du n {t,x)-Du{t,x),F)uv < {E^Du^xyDu^x)^) 1 ' 2 ^^^) 1 ' 2 -> 0, 

as n — >• oo. For the second term on the left-hand side of (|8"Tj) . we assume that 
F is smooth, i.e. i* 1 = F^ip with Fo € S and y; G W. By the Cauchy-Schwartz 
inequality and the Dominated Convergence theorem, 



E{G{t-;X-*){u n -i-u),F)- H 'p = aH I E[F (u n -i(s,y)-u(s,y))} 

J{R + xR d ) 2 

G(t -s,x- yMs', y r )\s - s'\ 2H - 2 f{y - y^dydydsds 1 



<a H E{F 2 ) 1 ' 2 f {E\u n ^{ S ,y)-u{ S ,y))\ 2 ) l / 2 G{t- S ,x-y)v{s',y') 

J(R + xR d ) 2 

\s - s'\ 2H ~ 2 f(y - y')dydydsds', as n -> oo. 
This concludes the proof of ([50]) . □ 
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